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Open charm meson spectroscopy: 
Where to place the latest piece of the puzzle 

Pietro Colangelo and Fulvia De Fazio 
Istituto Nazionale di Fisica Nucleare, Sezione di Bari, Italy 

We discuss how to classify the cs meson Z)sj(3040) recently discovered by the BaBar Collabora- 
tion. We consider four possible assignments, together with signatures useful to distinguish among 
them. 

PACS numbers: 13.25.Ft, 12.39.Fo, 12.39.Hg 
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In the infinite heavy quark mass limit, heavy-light Qq 
mesons can be conveniently classified in doublets labeled 
by the value of the total angular momentum = s g -I- £ 
of the light degrees of freedom (light antiquark and glu- 
ons) with respect to the heavy quark Q is the 

spin of q and £ the orbital angular momentum. For 
H. = (s-wave states in the constituent quark model). 



the doublet has S( 



1 



and consists of two states (de- 



noted as (P, P*), with P a generic heavy meson) having 
spin-parity jj^ = (0^, l^)i/2- p-wave states, with £ = 1, 
form two doublets: {PS,P[) with jj^ = (0+,l+)i/2, and 
(Pi, PI) with JP ^ (l+,2+)3/2. For £ = 2, d-wave 
states give rise to other two doublets: (Pi*,P2) with 



(l-,2- 



J3/2 



and (P:^*,P3) with Ji; = (2-,3-)5/2 



The construction proceeds further, and has been suc- 
cessfully applied to classify beauty mesons, as well as 
mesons with a charmed quark, albeit in the latter case 
the finite heavy quark mass corrections can be important. 

In the case of cs mesons, both the lowest lying states, 
both the radial excitations (whose members of the var- 
ious doublets we denote by a tilde) can be described 
using this classification scheme. Among the known cs 
hadrons, the two mesons Ds{l9m) and i:);(2112) fill the 
doublet with jf^ ~ (0~,l~)i/2, while one can identify 
(P>;o(2317),I?;,i(2460)) and (I?si (2536), (2573)) with 
the doublets Jj^ = (0+, l+)i/2 and Jj^ ^ (l+,2+)3/2 
respectively (the discovery of 13*0(2317) and D'^^{2Am) 
is rather recent 0, and their classification has been the 
subject of discussions [1]). 

More recently, other candidates have been added to 
this list: Dsj(2860), discovered by the BaBar Collabo- 
ration and Z^s J (2710), discovered by the Belle and 
BaBar Collaborations [3, Q, both observed in the DK 
final state. In the case of Dsj{21\Q), it has been possi- 
ble to determine the spin-parity: ~ 1^, studying its 
production in B decays. As for Dsj{2?>QQ), the assign- 
ments = 3~, with radial quantum number ti = 1 [6|, 
and = 0+, with n = 2 (i.e. the first radial excita- 
tion of D^j (2317)) have been proposed soon after the 
discovery. 

The latest piece of information comes from the BaBar 
Collaboration Q, thanks to a new analysis of DK and 
D*K final states. Dsj{271Q) and Dsj{2860) are seen de- 
caying to both DK and D*K, hence they have natural 
parity 1", 2+, 3~, • • •. The decay into D*K ex- 



cludes the assignment = 0+ for Dsj{2860). Other 
information comes from the measurement of the ratios 



BR{D, j{2710) D*K) 

BR{Dsj{27lO) DK) 

BR{D,j(2860) D*K) 

BR{Dsj{2860) DK) 



0.91 ±0.13,tat± 0.12,2^,4 



1.10±0.15,tat±0.19,j,3t 



(1) 



where D^*^K is the sum over the final states D'^*'>'^K'^ 
and D'^*^^ K% Q. Comparing these data with the pre- 
dictions in [3], it seems very likely that i3sj(2710) is 
the first radial excitation of P'*(2112), while the case 
of P'sj(2860) requires further study since the measured 
ratio in ([T|) is larger than the theoretical prediction 
BR{Ds.,[28m)-,D*K) ^ 
BR{D,j{28m)^ DK) ' "■^'^ 

In the same analysis, the BaBar Collaboration ob- 
served another broad structure in the D* K distribution, 
Dsj{30^0), with [3 



A/(i?,j(3040)) 
r(i?,j(3040)) 



30AA±8statitf)syst McV 

239 ± 35 statittt) syst MeV 



(2) 



Studies of angular distributions for this state have not 
been attempted, due to the limited statistics. Here we 
discuss possible classifications for this new cs meson using 
the only information available, the measured mass, width 
and decay mode, together with the full set of information 
concerning the other levels/doublets. 

Since £'sj(3040) decays to D*K and not to DK , it 
has unnatural parity, with possible = 1+, 2^, 3+, • • •. 
The lightest not yet observed states with such quantum 
numbers are the the two = 2~ states belonging to the 

3 5 

£ = 2 doublets, with — 2 ^^'^ ~ 2' denote them 

as Ds2 and D'*2, respectively (the case ~ 3+ corre- 
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sponds to a doublet with S£ = — , which is expected to be 
heavier). In the case of radial excitations, the identifica- 
tion with the states with n = 2, ~ 1+, and S( = — 



(the meson D'^i) or si ~ - (the meson Dgi) is possible. 

These four assignments are indicated in Table U which 
reports the classification of all the other known cs mesons 
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TABLE I: Known cs mesons organized according to sf and J^; the measured masses are indicated. Allowed possibilities for 
L>sj(3040) are displayed as Dsj(3040)?. 
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£is(1965) (0-) 


Dsj{2317) (0+) 


D,i(2536) (1+) 




D,j(3040)?(2-) 


J'' = + i 


0:(2112) 


I>sj(2460) (1+) 


I>:2(2573) (2+) 


Dsj(3040)?(2-) 


i3,j(2860)?(3-) 


(n = 2) 














(0-) 


(0+) 


D,j(3040)?(l+) 




(2-) 


J"'^ = sf + "1 


D,j(2710) (1-) 


Dsj(3040)?(l+) 


(2+) 


(2-) 


(3-) 



(Dsj(2860) has been assigned to the — 2 doublet, with 

a question mark since confirmation is needed) . 

Some indications about the masses of these states 
come from potential model calculations. For example, in 
Ref. [ll[ the spectrum of heavy-light mesons is computed 
in the framework of a relativistic quark model (RQM), 
with results: 



{RQM) _ 



3114 MeV 
3165 MeV 
2953 MeV 
2900 MeV 



(3) 



Notice that, if the identification of ZJ^j (2860) as the 
= 3~y2 state were experimentally confirmed, this 
would disfavor the assignment of Dsj{30A0) to its spin 
partner D*2 with Jj^ = 2~^^ , since a mass inversion in 
a spin doublet seems unlikely. For a similar reason, one 
would also disfavor the identification of £>sj(3040) with 
Ds2, although in that case the two mesons would belong 
to different doublets. 

The four classifications for Dgj (3040) in Table U can 
be discussed computing the allowed strong decays. To 
this purpose, we work in the heavy quark limit in which 
the various spin doublets are described by effective fields: 
Ha for sf = ^ (a = u, d, s is a light fiavour index); Sa 

and Ta for sf = 5^ and sf = respectively; Xa 
and X'a for the doublets corresponding to orbital angular 



momentum 

Ha = 



2, i.e. si 



and 




(4) 



7-)*U.I>' 

^2a 757.. 



Plau\j 2 



XT 



i + f 



7o- 



r,*' aj3 

P2a \/375 



with the various operators annihilating mesons of four- 
velocity V which is conserved in strong interaction pro- 
cesses (the heavy field operators contain a factor ^/rnp 
and have dimension 3/2). 

Let us consider decays with the emission of a light 
pseudoscalar meson. The octet of light pseudoscalar 
mesons is introduced considering the fields: ^ = e , 
S = f ^, and the matrix Ai containing tt, K and rj fields 
{U = 132 MeV): 



M = 



( 



V 



2 " 



K 



K+ \ 



2 " 



(5) 



■q J 



At the leading order in the heavy quark mass and 
light meson momentum expansion, the decays F — > HM 
{F — H, S, T, X, X' and M a light pseudoscalar meson) 
are described by the Lagrangian interaction terms [l^ : 



Ch 
Ct 

Cx' 



g Tr 
hTr 



HaSblt.lbA'ia 



h.c. 



h' 
Ax 



A. 



HaT^{iD^,A+^PA^,)bal5 
HaX^{lD^^+lpA^)bal5 

HaX'^^''[ki{D^,D,}Ax 



hiD^DxA, + D,DxAf,)]^yj5 



+ h.c. 
+ h.c. 

h.c. 



(6) 
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where Afj^ba = H^^'^m? ~ ^^M^^^a- is the chiral 

symmetry-breaking scale: we use = 1 GeV. Cs and 
Ct describe transitions of positive parity heavy mesons 
with the emission of hght pseudoscalar mesons in s— and 
d— wave, respectively, and g, h and h' are the effective 
coupling constants. On the other hand, Cx and Cx' de- 
scribe the transitions of higher mass mesons of negative 
parity with the emission of light pseudoscalar mesons in 
p— and /— wave, with coupling constants fc', ki and fc2. 

At the same order in the expansion in the light me- 
son momentum, the structure of the Lagrangian terms 
for radial excitations of the H, S and T doublets does 
not change, being only dictated by the spin-flavour and 
chiral symmetries, but the coupling constants g, h and h' 
must be replaced by g, h and h' . The advantage of this 
formulation is that meson transitions into final states ob- 
tained by SU{3) and heavy quark spin rotations can be 
related in a straightforward way. 

With the effective Lagrangians in Eq.© we can evalu- 
ate the strong decays of Dgj {3040) to a charmed meson 
and a light pseudoscalar one in correspondence to the 
four classifications in Table HI In particular, the ratio 



Ri = 



r(i^,j(3040) ^ D*Tj) 
r(D,,/(3040) ^ D*K) 



(7) 



(£)(*) A' = D^*'>°K+ + i:>(*)^A'0) can be computed, with 
results for the various assignments: 



RiiDsi) 

Rl{Ds2) 

RiiD%) 



0.34 
0.20 
0.245 
0.143 



(8) 



It is important to notice that the dependence on the ef- 
fective couplings cancels in the ratio. The spread among 
the various predictions permits to discriminate among 
the assignments, in particular between D'^i and D*!^. 

The mass of -Dsj(3040) is large enough to allow other 
decay modes. Decays to the members of the doublets 
with sf = OT sf = and a pseudoscalar meson are 
possible: these are the modes {Dq,D[)K and {Di,D2)K. 
As for channels with the r] in the final state, the only kine- 
matically allowed mode is UJo?? (for D'^iT] the available 
phase space is tiny). The features of such decay modes 
are different for the four considered assignments. At the 
leading order in expansion, the states D'^^ and Dgi 
can decay to D^K, 0*^7], D[K, DiK and D^K in p— 
wave; Ds2 and D*2 both decay to D^K, 0*^7], D[K in 
d— wave; D*2 can decay to DiK and in d— wave, 

while in the case of Ds2 the decay to D2K proceeds in 

s— wave and the decay into DiK is allowed at O (^:;^) 
in d— wave. 

Other kinematically allowed modes of Dsj (3040) are 
those with the emission of a light vector meson, specif- 
ically decays into DK* or Dgcj) [3- It is possible to 
describe such decay modes using an approach based on 



effective Lagrangian terms analogous to the one followed 
above. In the case of light vector mesons in the final 
state, the method has been developed in Ref.[l^ on the 
basis on the hidden gauge symmetry idea We refer 
to the review [isj for a detailed description; here we only 
mention that in this approach the octet of light vector 

mesons is described by the field = where 

v2 

is a hermitian 3x3 matrix analogous to ([5]), containing 
the light vector meson fields p^'", K*^, K*^, ws- 

The constant gy can be fixed to the value gy = 5.8 by 
the Kawarabayashi-Sukuzi, Riazuddin-Fayyazuddin rela- 
tions [H]. 

In Ref.[3l the effective Lagrangian terms describing 
the transitions S — )- HV and T — >■ HV are reported, S 
and H representing the heavy doublets Sa and Ha in Eq. 
(|^ and V denoting a generic light vector meson. They 
read: 



SHV 



L 



THV 



zC5Tr[5,i/fc7^(V^-p^)ba 



(9) 



iprTr HaT^v" Ffj,^{p) 



ba 



where V^ba = f (C^^m^ + ?^A'?^)(,a' ^^'^ '^s,t, Ps,t are 
effective coupling constants, and the field strength tensor 
F^^ is defined as Ff_,^{p) = d^p^ - d^p^ + [p^, p^]. 

The calculation of all possible decay modes for each one 
of the four classifications of D^, 7 (3040) would require the 
values of several coupling constants, which are unknown. 
Nevertheless, some conclusions can be drawn. 

One can estimate the widths of the decays into DK* 
and Ds(j) for the two = 1+ states by the effective 
Lagrangians in PT7)1 using the values of ( and p: Cs = 
0.10 and ps = —0.10 GcV^^ computed for the analogous 
transitions involving the heavy mesons belonging to the 
lowest-lying doublet p^ . This provides us with hints 
about the contribution of such modes to the total width 
of Z^s J (3040) in the case it is either D'^-^ or Dgi (in the 
case of Ds2 and 0*2 these decays proceed in p-wave and 
hence they are expected to contribute less significantly): 



T{D',, ^ DK*) 

T{D,, ^ DK*) 
T{D,i ^ A,0) 



95 MeV 
44 MeV 
15 MeV 
6 MeV . 



(10) 



It can be noticed that the contribution of these modes to 
the full width should be sizable only in the case of D'^^^; 
therefore, this assignment can be disentangled through 
the experimental analysis of the DK* and Ds(j) decay 
channels. 

After having discussed four possible classifications for 
Z)sj(3040), it is interesting to consider some features of 
its spin partner in the various cases [20| . 
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TABLE II: Features of the decay modes of _Dsj(3040) and of its spin partner at leading order in expansion for the four 
proposed assignments. 



decay modes 


Di, {n = 2, Ji^ = 1+2) 


Dsi {n = 2, Ji^ = 1+2) 


Ds2 (n=l, Ji^ =23-3) 


(n = i, Ji^ = 2-/2) 


D*K, Din 


s— wave 


d— wave 


p— wave 


/— wave 


Ri 


0.34 


0.20 


0.245 


0.143 


Dr,K, D^oV, D[K 


p— wave 


p— wave 


d— wave 


d— wave 


DiK 


p— wave 


p— wave 




d— wave 


DIK 


p— wave 


p— wave 


s— wave 


d— wave 


DK*, Ds4> 


s— wave 


s— wave 


p— wave 


p— wave 


r ~ 140 MoV 


r ~ 20 MeV 


neghgible 


neghgible 




spin partner 




Dlo {n = 2, Ji^ = 0+2) 


£>:2 (n = 2, Ji^ = 2+2) 


D'si {n = 1, Ji^ = 13^2) 


Ds-s {n = 1, = 35/2) 


DK, DsTj 


s— wave 


d— wave 


p— wave 


/— wave 


D'K, Dlri 




d— wave 


p— wave 


/— wave 


d:,k, Dion 






d— wave 




D[K 


p— wave 


p— wave 


d— wave 


d— wave 


DiK 


p— wave 


p— wave 


s— wave 


d— wave 


d;k 




p— wave 




d— wave 



• If i:»s,7(3040) is D'^i (sf = = 1+, n = 2), its 
spin partner is D*q, a J-^ = 0+ state, the first radial 
excitation of Dsj{2317). This state can decay to 
DK and DgT] in s-wave; p-wave decays to D[K 
and DiK are also allowed. 

• If D,j{3040) is D,i {sf = | + , = 1+, n = 2), 
its spin partner is D*2 with = 2+. It is allowed 
to decay to DK, DsTj, D*K, D*r], DK* and Dscj) 
in d-wave, and to D[K, DiK and D2K in p-wavc. 

• If L>^j(3040) is D,2 (sf = f, = 2", n = 1), its 
spin partner is the vector meson D*i with = 1^. 
It can decay to DK, DsT], D*K, D*'q, DK* and 
Dgcj) in p- wave, to DqK, -D*q?7 and D[K in d— wave 
and to DiK in s— wave. The decay to D2K is 

allowed at O ( — ) in d— wave. 

• If £'.7(3040) is D*!^ (sf = = 2", n = 1), 
its spin partner is Dss with J-^ = 3^, decaying to 
DK, DsT], D*K and D^i], DK*, D^cf) m /- wave, 
and to D[K, DiK and D2K in wave. 

Since £'s,/(3040) has a broad width, we expect that also 
its spin partner shares the same feature. Considering the 
previous list, we can argue that D*q is broad due to its 
s-wave decays into DK and DsTJ. Also D*i has allowed 
s-wave decays, but only to DiK which is suppressed by 
phase space effects. 

The identification of -D.j (3040) with D'^^ is supported 
in Refs-fTrl. [Tsj on the basis of the cs mass spectrum 
[l7j or of the decay widths computed in the ^Pg model 



[l8|. In the second case, the identification with D'^i and 
Dsi is discussed: the full widths of these two states are 
computed and compared to the experimental measure- 
ment of r(_Ds,7(3040)), concluding that for D'^^ the ex- 
perimental width can be reproduced, with the predic- 
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FIG. 1: Spectrum of the cs system. All observed Dsj states, 
with mass indicated on the y axis, are assigned to a level 
with and proper name. The four assignments discussed 
for _Dsj(3040) are shown in correspondence to the mass value 
M = 3040 MeV. 
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tion: 



r(i:i^j(3040) -> DK* +Ds 



0.79. In our dis- 



r(i:i^j(3040) -> D*K + D*ri) 
cussion we have examined other information that can be 
exploited. 

To summarize the results of our considerations and 
compare the distinctive features of the four proposed as- 
signments for Ds J {30A0) (illustrated in FigH] in which 
all the known cs states have been included), we collect 
our findings and observations in Table [III In particular, 
we emphasize the role of the final states DK* and Dgcj), 
which deserve an experimental investigation. Moreover, 

3 ^ 

although the identification with a = 2^ = — 



state seems less probable, this assignment can be dis- 
carded/confirmed studying the s—wave final state. 
Search of the spin partner in each doublet would provide 
further information, enriching the cs spectrum recently 
disclosed by experiments. 
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